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NOMENCLATURE 


= Volume average crystallite diameter, nm 

Kq = Crystallite shape factor, dimensionless 

K,K' = Homogeneous reaction equilibrixm constants 

for the reactions, meta to ortho & meta to para 
isomers respectively, dimensionless 

A = meta xylene 

R = ortho xylene 

S = para xylene 

Rate constants for adsorption of A, R or S 
on catalyst surface, g,mol/(hr) (g. catalyst ) 

k,]c' = overall rate constants for the reactions, meta to 

ortho & meta to para respectively, g.mol/ 

(hr) (g. catalyst) 

Ka/Kr, K g= Adsorption equilibrium constants for components 
A, R,S respectively, (atra“l) 

= Rate of isomerization of m-xylene to p-xylene, 
g.mol/(hr. (g. catalyst) 

rg = Rate of isomerization of m-xylene to p-xylene, 

g.mo 1/ (hr ) (g- catalyst ) 

r^ = Rate of reaction of any particle size, g.mol/ 

(hr) (g. catalyst) 

r 2 = Rate of reaction of powdered catalyst, g.mol/ 

(hr) (g. catalyst) , 

-1 

LWHSV = Liquid weight hourly space velocity, hr 

P = m-xylene feed rate, g, mol/hr 

W = Weight of catalyst, g 

Specific metal surface area 


s 



X 


= Conversion of m-xylene to o-xylene 

= Conversion of m-xylene to p-xylene 

Greek Symbols t 

P 1, = Width of nickel pealc at half maximura intensity^ 

radians 

= VJavelength of X-ray radiation/ nrn 
% - Effectiveness factor/ dimensionless 

Q - Bragg ' s angle, radians 
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ABSTRACT 

Kinetics of m-xylene isomerization has been studied in 
vapour phase over nichel hydrogen mordenite (Ni^HM) catalyst. 
Catalyst was prepared by impregnation technique and tlae reactioi 
was carried out in hydrogen atmosphere in a fixed bed flow 
reactor. Niclcel content of the catalyst was varied from 0.5% 
to 5% by weight. Mass transfer effects were found to be negli-| 
gible. Activity and selectivity of the catalyst were affected ! 
by H 2 /m-xylene mole ratio, hquid weight hourly space velocity ; 
(LWHSV), reaction temperature, nickel content and nickel crys- ; 
tallite size. 1.5% nickel hydrogen mordenite gave maximum i 

conversion and activity. 

A lepara'lleiv: reaction scheme v/ith no interconversion 
between ortho & para isomers was, proposed for the reaction. 

The kinetic data were analysed on the basis of Langmuir~ 
Hinshelwood kinetics v;ith statistical data interpretation. 

The rate controlling step in the over all kinetics of the react* 
was dual site surface reaction model. The activation energy 
values wereI7.6lkcal/raol for p-xylone formation andl5.28 kcal/ 
mol for o-xylene formation respectively. . I 
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INTRO DUCTION 


At. the start of the 20th century the only use of xylene 
was as a component of solvent naphtha* At around 1920 the 
xylenes found their utilization in the dye industry as a compo- 
nent for lacquers. During this period the antiknock properties 
of xylenes were discovered. The usage of individual isomers 
ortho, meta and para xylenes were not recognized until the 
middle of 20th century when methods of separation of the indi- 
vidual isomers were developed. In 1955 the commercial utiliza-i 
tion of the isomer para xylene started as the raw material for J 
production of polyester fibres and films. In: 1960, the use of 
ortho xylene began in the manufacture of alkyd resins, an 
important component in making hard durable lacquers. The main ; 
component of these resins, phthalic anhydride could be obtained! 
by oxidation of naphthalene or ortho xylene. Now, vfith the 
great reduction in the quantity of naphthalene available, the I 
demand of ortho xylene for this purpose has increased many fold]. 
Another important use of phthalic anhydride is as an intermediai 
in the manufacture of plasticizers. 

The above mentioned uses along with many others placed 
the para and ortho xylenes as a basic raw material in a key 
position of importance in the petrochemical industry. 
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It is a paradox that m-xylene has got the least demand 
inspite of it being present in the largest amount as a product 
from the catalytic reforming of petroleum naphtha, Meta xylene: 
can be oxidized to isopthalic acid which is used for the 
production of unsaturated polyesters but the tonnage employed 
is very small. The bulk of the meta xylene has to be used as 
a blend in motor gasoline. So, as far as economics and demand ; 
of the isomers are concerned, if m~xylene can be isomerized 
successfully to o- and p~xylenes, it will be of great advantage 

The isomerization of xylenes is generally carried out | 
in the presence of an acid catalyst (1-3), The reaction can be| 
affected in both liquid and vapour phase. The liquid phase I 
isomerization using homogeneous acid-halide catalysts suffer 
from corrosion and recovery problems . The vapour phase isomeril 
zation on heterogeneous silica alumina and dual function type 
catalysts avoids these problems of homogeneous catalysts but 
the high temperature required leads to coke formation, thereby 
affecting the yield of the desired product by poisoning the 
catalyst. Addition of hydrogen to the xyl.;;ne feed reduces cokel 
formation, suppresses side reactions and increases the rate of 
xylene isomerization. 

In recent years, the aromatic isomerization reactions 1 
are studied using various zeolite catalysts because of their i 
better activity and selectivity compared to silica alumina and 
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homogeneous acid halide catalysts. Metals were dispersed on 
these supports to obtain better yield and selectivity of products 
Platinum has been extensively used for this purxjose. 

Hydrogen, mordenite has been reported as a catalyst for 
xylene isomerization in liquid phase (4,5) and vapour phase (6).: 
Nickel supported catalysts were used for isomerization of para- 
ffins (7), In view of these facts it is possible that nickel 
supported over mordenite may show considerable catalytic activit] 
for the isomerization of m-xylene. 

The aim of the present study vras to obtain xylene isomeri| 
zation rate data using nickel hydrogen mordenite (Ni-IiM) as the ; 
catalyst and to compare these data with reported values in the i 
literature for different catalysts. The isomerization reaction 
was carried out in the vapour phase using nickel supported 
mordenite catalyst. Nickel was supported on hydrogen mordenite 
by impregnation technique and reduced to nickel metal in the I 
presence of flowing hydrogen gas. The reaction was carried out 
in a fixed bed reactor under hydrogen atmosphere to keep the ; 
catalyst free from coke formation throughout the reaction. - 
Effect of liquid weight hourly space velocity (LWHSV), reaction I 
temperature, nickel concentration and nickel crystallite size | 
were the important variables studied for isomerization activity 
and selectivity of the catalyst. The kinetic data were analyzed 
using a model based on Langmuir Hinshelwood kinetics. 



CHAPTER 2 


LITERATURE SURVEY 

Xylene isomerization has been studied by a nuitiber of 
research workers in the past. Various types of catalysts v/ere 
employed for the reaction. These are briefly discussed below. 

The isomerization of xylene is greatly aided by the use 
of an acid catalyst and can be carried out in either the liquid: 
phase or vapour phase. Thermal isomerization of xylenes is of 
very little interest because extreme conditions are required and 
the yields are very poor (8). The liquid phase isomerization of 
xylenes can be carried out under either homogeneous or hetero- 
geneous conditions with a combination of Lewis acids as a 
catalyst/ for example, HP/BP^. Such reactions have been observec 
to be first order with respect to xylene concentration (2,3) 
and first order with respect to catalyst concentration for 
homogeneous phase conditions (2). 

Most of the earlier stiidies of liquid phase isomerization I 
of xylene were carried out using homogeneous acid halide catalysi 
M'claulay and Lien (1) studied the isomerization reaction using 
hydrogen bromide and boron trifluoride as catalysts and found I 
that the extent of isomerization reaction strongly depended upon 
the boron trifluoride concentration. Brown and Jungk (2) reports 
the values for rate constants and activation energies for the 




isomerization reaction using hydrogen bromide and aluminium 
bromiide as catalyst. They also proposed that the reaction 
would be occuring via series mechanism i,e., p~xylene 
m-xylene o*-xylene. The reaction followed second order 

kinetics, first order with respect to both xylene and catalyst 
concentration. Similar results were reported by Allen and 
Yats (3). The disproportionation of xylenes to toluenes and 
trimethyl benzenes, a major side reaction was suppressed using 
toluene as a diluent. The usage of homogeneous acid-halide 
catalysts in the liquid phase isomerization suffers from various 
drawbacks such as side reactions, sludge formation, loss of 
partially soluble catalyst, corrosion and difficulty in recovery 
of products. 

Due to the above problems of liquid phase catalyst- 
reaction system interest has turned towards acidic type hetero- 
geneous catalysts such as silica alumina. Unfortunately, these 
catalysts are not very active at the ordinary temperature and 
therefore high temperatures are needed for these catalysts to 
achieve considerable amount of catalytic activity. Most of the 
studies were carried out in vapour phase. The vapour phase 
isomerization of xylene over silica alumina catalysts were 
studied by number of workers (9-17), 

Hanson and Engel (16) reported the kinetics of vapour 
phase isomerization over silica alumina catalyst. They found 
that m-xylene isomerizes reversibly to o-xylene and p— xylene 
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over silica aliomina/ but direct interconversion between o- 
xylene and p~xylene does not occur. Xylene isomerization was 
found to be first order and was described by a single site rate 
model. The values of rate constants and activation energies were 
reported for isomerization as well as disproportionation. The 
activity of silica alumina catalyst for isomerization of xylenes ^ 
was found to be strongly dependent on the combined effect of 
time/ temperature and gaseous environment. Cortes and Corma 
(17) investigated the mechanism and kinetics of vapour phase 
isomerization of xylenes over silica aluraina catalyst. Silvestr: 
and Prater (15) investigated the kinetics of the selectivity 
of gas phase xylene isomerization over silica alumina catalyst. ; 
They proposed a coupled reaction scheme involving six rate 
constants 


^ p-xylene 
o-xylene 

where they considered the interconver s ion among the three 
isomers. The coupled reaction network is difficult to solve 
since it is seldom possible to obtain absolute rates for the 
three reactants due to catalyst aging/ unknown catalyst site 
density etc. Fortunately/ a method developed by Wei and Prater 
(18) permits reduction of the coupled reaction network to an 
uncoupled network. This mathematical method permits one to 




m-xylene 
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calculate the six rate constants without the need to obtain 
absolute conversion rates for any of th o three reactions. 

Silica alumina alone is not very efficient for the _ 
isomerization reaction* To enhance the activity and selectivity 
of the reaction, various metals supported on silica-alumina are 
used as catalysts . Platinum is extensively used for this 
purpose. Isomerization activity of the xylenes over a platinum 
containing silica alumina catalyst was improved by coupling 
with a recycle of product system (13). Isomerization of o- and 
p~xylenes over a bifunctional Pt/Si 02 -Al 202 catalyst studied 
in a temperature range of 703K to 748K indicated that m-xylene 
was an intermediate in both the forward and backvrard reaction 
(11). A catalyst containing platinum, rhodium and tin over 
iT-alumina showed more than 80% conversion of m-xylone to 
p-xylene ( 14 ) ,. 

Activities of various zeolite catalysts were established 
long back for the paraffin isomerization reaction. These 
catalysts have shown higher activity and selectivity for the 
aromatic isomerization reaction. Infact, in recent years, most 
of the xylene isomerization studies have been carried out on 
various zeolite catalysts which showed considerable improvements 
over the silica alumina and the homogeneous acidic halide 
catalysts as far as activity and selectivity are concerned. 

Synthetic zeolites or molecular sieves have a general 
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foarmula as % 

2^3 * * ^^2 ^ 

where ^ M is a cation of valence n. The framevrork structure 
contains channels and Interconnected voids which are occupied 
by the cation and water molecules. The cations are quite 
mobile and may be exchanged to varying degrees by other cations 
The intercrystalline water can be removed reversibly : for some 
zeolites. The mordenite framework structure is shown in Pig. 1 
There are various types of zeolites but the most commercially 
useful zeolites are A type, X type, Y type and mordenite. 

The unit cell representation of these zeolites are (35) s 


Type 

A s 

Nai2 

[(A102)i2 

(SiO^ 

^12 J • 

27 

H 2 O 

Type 

X s 


[_(A102)86 

(Si02 


264 

H 2 O 

Type 

y s 

Nase 

[_(A102)56 

(Si02 

^ise]* 

250 

H 2 O 

Mordenite 

sNag 

1 (A102)8 

(Si02 

^40 1 * 

24 

H 2 O 


The molecular sieves can be located on the three phase 
(silica-alimnina-alkali ) compositional diagram as shown in 
figure 2. in general, the acid resistance increases as the 
silica content increases. Molecular sieves based on mordenite 
material with a silica-alxxmina ratio of 10 s 1 have much more 
chemical and thermal stabilities than those having lower silica 
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1007, SiO,.' 



FJG.2 GENERAL CHEMICAL COMPOSITION OF 
ZEOLON MATERIALS 



-IG-I MORDENITE FRAMEWORK STRUCTURE 

eo -PLANAR FACES 
L-] UNIT 





10 


aliimina ratio. 

A number of studies are reported in the literature for 
isomerization of xylenes over zeolite catalysts. Hansford and 
Ward (19) found higher isomerization activity of zeolite cata- 
lysts compared to silica-alumina catalysts. Catalyst acidity 
played an important role in the activity of the catalyst. 
Spectroscopic studies indicated that the catalytic activity of 
aluminosilicates for m-xylene isomerization depends on the 
concentration of strongly acidic centres (20). Bremar et al. (21: 
revealed that the catalytically active centres in xylene isome- i 
rization are Bronsted acid centres which occur on cation exchan-t 
ged zeolites by dissociation adsorption of H 2 O molecules if the 1 
polarizing effect of the cations is sufficiently strong and : 
the cations are not saturated sterically or bonded by the 
zeolite skeleton. Experiments conducted by Sharma et al. (22) 
on o-xylene isomerization with various forms of 2 c-zeolites j 

showed that calcium-X and hydrogen-X possessing stronger acid 
sites were more active than sodixm-X which contains only weak I 

i 

acid sites. 

To obtain greater activity, the acidity of the zeolites | 
can be increased by increasing the silica/alumina ratio. This ] 
can be achieved by dealuminating the zeolite. Gajewskiet al. (2 
obtained higher activity of the Al-deficient zeolite (dealumina 
Y-zeolite) in m-xylene isomerization. Parasiewicz et al. (24) j 
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used series of X and Y zeolites in the Lanthanum (La), 
calcium (Ca), magnesium (Mg) and hydrogen (H) forms as catalyst 
for mc-xylene isomerization. Best activity was exhibited by 
dealumina ted mordenite containing 20 si silica-alumina ratio 
in the hydrogen form which gave around 40% para and ortho 
xylenes with 80% selectivity at 473 °K. The pore structure 
of mordenite was found to be an important factor for the 
process studied. Papp and Miklosy (25) found that the cata- 
lytically active centres of hydrogen mordenite .for xylene 
isomerization reaction were the OH groups. 1/2 methyl group 
shift was found to be the mechanism of reaction. Hopper. and 
Shigemura (4) obtained good isomerization activity and selec- 
tivity for o-xylene isom.erization in liquid phase over hydrogen 
mordenite catalyst using toluene as diluent. First order 
kinetics was applied to fit the model satisfactorily. Norman 
et al. (5) investigated the kinetics of liquid phase 

isomerization of xylene over hydrogen mordenite, A Langmuir- 
Hinshelwood kinetic model was developed for the isomerization 
of xylene wliich assumes first order interconversion among the, 
three isomers. This model is very comprehensive as it includes! 
the effects of feed composition, catalyst activity, and 
temperature. The effect of composition was included in the J 
reaction rate, equation by the Langmuir-Hinshelwood expression ; 
which contains separate adsorption parameters and reaction rate 
parameters. The value of activation energy was reported. : 
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Lanewala & Bolton (26) proposed that xylene isomerization occured 
by the transalkylation mechanism. The isomerization of xylenes 
over a zeolite catalyst is invariably accompanied by transalky- 
lation and the degree of isomerization is proportional to the 
extent of transalkylation. They proposed that the reaction 
proceeds via a bimolecular mechanism involving diphenyl type 
transition state. Using the monomolecular rate theory developed 
by Wei and Prater , Chutoransky & Dvryer (27) analyzed the kinetics 
of the liquid phase isomerization of xylene over a zeolitic 
catalyst. The kinetic analysis was presented primarily in 
term.s of the time independent selectivity kinetics. It was 
found that the isomerization is a simple series reaction | 

o-xylene m-xylene p-xylene. Collins et al. (28), i 

investigated xylene isomerization and disproportiontion over [ 

Lanthanum Y catalyst. A coupled reaction scheme involving six 
rate constants was proposed. The relative rates of the isomeri- ; 
zation reactions were calculated using Wei-Prater method. The 
results were found to be consistent with a series reaction path 
for the isomerization reaction. Collins and Medina (29) studied j 
xylene isomerization using ZSM-S zeolite catalyst. They found J 
that disporportionation was not a significant reaction in the I 
tenaperature range of 523-623K, The relative rate constants 
calculated using VJei-Prater technique appeared to be influenced 
by diffusion and/or shape selectivity. ZSM-5 catalyst showed 
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a slight selectivity for p-xylene formation during xylene 
isomerization. This was believed to foe duo to shape selective 
diffusion. 

i'iiklosy et al. (6) studied xylene isomerization on H- 
forrns of natural mordenite and clinoptilo.lite, prepared foy one 
of the following procedures 2 (i) acid treatment (ii) NH 4 '^ 
exchange followed foy daamrnoniation and (iii) Ag ' exchange 
followed by reduction with hydrogen. They reported that the 
active sites for xylene isomerization are the acidic OH groups 
on the outer crystal surface of the H-forms of mordenites and 
clinoptilolite prepared differently. Bronsted centres of 
different activity/ selectivity and stability are formed upon 
treatment with hydx'ochloric acid, or ammonium ion-exchange 
followed foy thermal deairimoniation/ or silver ion exchange follo- 
wed by reduction with hydrogen, H-forms of different zeolites 
prepared similarly showed different catalytic behaviour/ thus 
indicating the effect of the lattice structure. The highest 
activity without deactivation was attained vrith the H(Ag)-form 
of natural small part mordenite/ which contains accessible 
Bronsted acid sites in maximum concentration with maximura 
efficiency and stability. Thus the kinetics of isoraerization 
could be investigated reliably on this catalyst. 

On this catalyst the isomerization is a consecutive 
monomolecular transforrnatiian corresponding to the scheme 
0“xylene - - — ^ ra- xylene p-xylene. The kinetic order 
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varies between zero and unity, indicating a medium coverage of 
active sites. In terms of the Langmuir-Hinshelviood Jcinetics, 
the rate determining steps are the surface reactions, in which 
the shift of methyl groups to neighbouring positions take 
place in the adsorbed species, which is presiomably a benzenium 
ion. The reactivities of xylene isomers do not differ substan- 
tially. The adsorption equilibrixom constants of the isomers are 
in the same order of magnitude and their values affect the 
conversion rates only to a small extent. 

Recently much of research works on xylene isomerization 
is directed towards getting high para selectivity feat this 
reaction. Ibstly investigated catalysts in this field are the 
ZSM-5 type. These catalyst supports have pore dimensions 
greater than the critical 'dimension of p-xylene but very close | 
to those of m- and o-xylenes thereby lowering the m- and o- | 

xylene diffusivities to a value much lower than that of p-xylenoJ 
and thus giving rise to high para selectivity, A ZSM-5 catalyst I 

(29) was found to yield a para/ortho ratio of 1.6 in the iso- 

merization of m-xylene. A significant extent of disproportionat: 
was also observed. A HZSM-5 catalyst containing platinum & I 

Tin (30) gave para selective isomerization of mixed' xylenes . I 

But at the same time some amount of hydrocarbons, benzene ■ 

and toluene were formed as side products. The tin present in 
the catalyst was believed to suppress hydrogenation of the 
benzene ring, i 
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Use of nickel raetal loaded zeolites for the isomerization 
of n-paraffins have been studied for a long time. Burch (7) 
studied the isomerization of n-hoxane on Ni/Al 202 catalyst. 
Minachev et al, (31) reported good activity for isomerization 
of butanes over nickel-hydrogen mordenite catalyst. Beecher 
& Voorhies (32) reported paraffin isomerization over a hydrogen 
mordenite based catalyst and concluded that the mordenite based 
catalysts are more active than faujasite based catalysts for 
these reactions. Similar observations were made by Inou & 

Sato (33) for xylene isomerization reaction. They found 
better activity and selectivity of mordenite based catalysts 
than amorphous silica-alumina & Y-type zeolites for isomeri- 
zation of xylene. 

Sonsarma (34) studied the isomerization of m- xylene 
in the liquid phase over nickel hydrogen mordenite, 1.5% nickel 
hydrogen mordenite gave the maximum conversion and activity. 

It appears from the literature that nickel hydrogen 
mordenite catalyst has not been studied extensively for vapour 
phase m-xylene isomerization. It is believed' tha.t this catalyst 
will show better catalytic activity and higher isomerization sole 
tivity. Moreover, limited information is available regarding the 
kinetic study of this reaction using mordenite as a support. A 
detailed study of m-xylene isomerization including its kinetics 
has been undertaken using nickel supported mordenite i 
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C HAPTER 3 
EXPERII-4ENTAL 

^ ^ Catalyst supports and chemicals 

Zeolite hydrogen mordenite (Zeolon 900 H) was obtained ; 

from M/s. Norton Ci'iemicals, England. Typical physical characterise, 
of this zeolite are t 

10 I 

H 

8-9 A I 

0.64-0.72 g/cc 
2 

400 - 450 m /g 

Pure grade m-xylene was obtained from M/s. Reidel-Dc HaSNAG Seelzel 
Hannover . The pPrity _ of m-xylene was checked by GLC 

and was more than 99.5%. Anal R grade nickel nitrate was 
obtained from hl/s, Thomas Baker & Co., Bombay. Technical j 

grade hydrogen and nitrogen (99.5%) gases were supplied by 
^i/s. Indian Oxygen Limited, Kanpur. I 

3.2 Cataly^st prepara^^ i 

The nickel-hydrogen mordenite (Ni-HM) catalysts 
containing nickel metal on hydrogen mordenite v/as prepared by | 

impregnation technique (36) as described below, f 


Si02/Al202 % 

. Cation s 

Effective pore diameters 
Bulk packing density s 
Surface area s 
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The amount of vrater absorbed by the support material 
(grams water absorbed pei' gram support) was first determined. 
Standard nickel nitrate solution was prepared using double 
distilled water. The amount of standard nickel nitrate solution 
to be impregnated was calculated for the preparation of a 
catalyst of known nickel percentage. This total amount of 
solution was impregnated in several instalments . In each 
instalment/ appropriate amount of solution, equal to the 
amount that can be absorbed by the support was added in a 
rotary evaporator. The rotary evaporator was maintained at 
a constant temperature of 383 t IK. After the requisite 
amount of solution v/as added the material was dried in an 
oven maintained at 393 + 5K. The catalysts were subsequently 
calcined at 823 K in a muffle fiarnace for 20 hrs. The 
catalysts were rediiced in situ in flowing; hydrogen (40 ml/min) 
at a temperature of 673 K for 15 hrs. The reduced catalysts 


were cooled to room temperature in hydrogen atmosphere. 


3 . 3 Experimental 


Isomerization of m-xylene was carried out in a fixed 
bed flow reactor. A schematic diagram of the setup is sho’t'jn 
in Figure 3. The hydrogen/m-xylene mixture, used as reactants, 
was obtained by bubbling hydrogen through two bubblers, 
arranged in series containing m-xylene kept in a constant 
temperature bath at a temperature of 364 + IK to have a 
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desired hydrogen/ra-xylene ratio. The reaction mixture was 
preheated in a tubular furnace which was controlled v;ithin 
an accuracy of + 5K by an IBP temperature controller using 
iron constantan ' thermocouple. The reaction mixture 

was passed through the reactor. The reactor was heated in a 
tubular furnace. The temperature of the catalyst bed was 
controlled within an accuracy of + IK by an APLAB temperature 
controller using chromel-Alxoinel thermocouple placed centrally 
in the reactor. I 

3,4 Expeirimental procedure 

3.4.1 Activity and selectivity study s | 

I' 

3.0 g of a catalyst saaiple was charged in the reactor, | 
The catalyst was reduced in situ in the reactor in flowing I 

hydrogen (40 ml/min) for 15 hrs, before any activity measxare- I 

ments were made. The reaction .temperature was adjusted to the 

' ‘ , 'I 

desired value. The hydrogen flow rate v/as set at 32 ml/min I 

and allowed to stabilise for naif an hour. The hydrogen/m- 
xylene mixture was allowed to pass through the catalyst bed. 

Liquid samples were collected after 2 hrs once the steady state I 
was attained, 

3.4.2 Kinetic study s 

i' 

Kinetic study was performed under differential | 

conditions using 1.0 g of activated catalyst sample. The | 
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hydrogen flow rate used was 43 ml/min and allowed to stabilize 
for half an hour. The kinetic data were taken at different 
temperatures ranging from 523-623 K under different partial 
pressures of -m-xylene. The partial pressiore of m--xylene was 
changed by introducing nitrogen as a diluent. Nitrogen was 
mixed with the hydrogen/m-xylene mixture and allow’ed to pass 
through the catalyst bed. Liquid samples were drawn at interva 
of 2 hrs and analyzed in GLC for product composition. 

3.5 Prod uct a nalysis . 

Analysis of the liquid samples drawn v/ere performed 
using a CIC Gas chromatograph fitted with a flame ionization 
detector. The components of the sample mixture were separated 
over a 1.6 mm i.d. stainless steel colmin having a length of 
3 m and packed with 5% Bentone - 34, 5 % Di-isodecyl pthalate 
on 6 0-80 mesh cliromosorb W, Pure nitrogen v;as used as carrier 
gas. Carrier gas flow rate of 15 ml/rain, hydrogen (fuel) gas 
flow rate of 15 ml/min and a column temperature of 383 K were 
found to be the best operating conditions for optimum separation 
of products. A typical chromatogram showing peaks of toluene, 
p-xylene, m-xylene, o-xylene and trimethylbenzenes is shown 
in Fig. 4 as a function of retention time. Mixtures containing 
known amounts of the above compounds were analyzed under the 
same column conditions. The areas under the peak of the 
components were found to be directly proportional to the 




FIG.4 TYPICAL CHROMATOGRAM SHOWING PRODUCT DISTRIBUTION FOR META XYLENE 

OVER 1 S^/o Ni-HM CATALYST AT 598 K 



22 


respective mole % of the components. Therefore the composition 
of any component was equated to the ratio of the area under the 
peak of that component to the Siam of the areas of all the peaks 
appeared. 

^ ^ ?r.^?Z...i3-if fraction studie^s 


The X“ray diffraction pattern of various nickel catalysts 
and the untreated support material was obtained on a horizontal 
type Siefert diffractometer, employing CioK.g^ radiation. This 
was done in order to check the loss of crystallinity of the 
catalysts during pretreatments and reaction conditions. The 
average nickel crystallite size was also determined for 
catalysts of various nickel concentrations by the X-ray line 
broadening technique (37). The X-ray diffraction data for 
the catalyst samples and the zeolite suppox't were recorded 
over the range of 20°“50° at a goniometer speed of 1.2'^ /min 
in 20 (0 = Bragg'fe angle), a chart speed of 15 mm/min and 
counts per second of IK. The nickel (ill) peak appeared at 
44.52'= in 20. 

The average nickel crystallite size was determined 
using the relation given by Klugh and Alexander (37 )i 


V 


l&i^ Gos 0 


( 1 ) 


where, d^ = 


average crystallite diameter, nm 
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Kq = crystallite shape factor = 0,9 

71 = wave length of CioK.,-. radiation = 0« 154178 nin 

R 1 = width of nickel peak at half maximum 
r h 

intensity, radians 

0 = Bragg's angle at which the Ni(lll) peak appeared 

= 22,26°' with C'uK^_ radiation 


The specific metal surface area was also calculated from XRD 
measurements with the assT.imption of spherical metal crystallites 
using the relation given below 5 


s 


■where, 0 ~ 


d^ = 


10^ mVg 

J V 

bulk density of nickel metal - 8.90 g/ml 
average crystallite diameter, nm 
specific metal surface area, m^/g 


( 2 ) 


s 
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CHAPTER 4 

RESULTS AInTD DISCUSSION 


Experimental runs were performed over a series of nickel 
hydrogen mordenite (Ni-HM) catalysts in vapour phase using pure 
m-xylene in a fixed bed flow reactor. The present study was 
conducted in two parts s 


(a) Activity and selectivity, studies 

(b) Reaction kinetics and modelling 

Catalytic conversion, activity and selectivity are 
defined as follows s 


Conversion 
(mole %} 


moles of m-xylene reacted 

moles of m-xylene fed in the reactor 


X 100 


(3) 


Isomerization 
activity (mole %) 


moles of m-xylene isomerized i qq 

moles of m-xylene fed in the reactor 


(4) 


Isomerization 
selectivity 
(mole %) 


moles of m-xylene isomerized 
moles of m-xylene reacted 


X 100 


(5) 


Since p-xylene is an iirportant product, the following 


terms are also defined s 
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p-xylene 
selectivity 
(mole %) 


moles of p-xylene formed ^ 
moles of m- xylene isomeri'sed 


p-xylene yield = 

(mole 3%) moles of m-xylene fed initially 


(7) 


Experimental runs were also taken with blank hydrogen 
mordenite support under similar reaction conditions to check its 
catalytic activity. The reproducibility of experimental data 
was checked by repeating the experiments under exactly same 
reaction conditions. The results were found to be reproducible 
within an error of + 5%. 

In a kinetic study it is very important that mass transfer 
effects be absent during the reaction in a fixed bed catalytic 
reactor. Mass transfer effects, both external and pore diffusion 
were checked in the present study before activity, . selectivity 
and reaction kinetics were studied. 

4.1 Mass Transfer Effects 

The important variables which can affect the kinetics 
of the reaction due to diffusion resistance are s 

(a) Space velocity (LWHSV) . 

(b) Catalyst particle size 
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4.1.1 Effect of space velocity (LVjHSV) 

Effect of external diffusion on conversion was studied 

T , . W (q-. catalyst) (hr) 

by conducting two runs at varying ^ ■= — — (gniioTl 

with constant weight of catalyst. The weight of catalyst in 

one run was double the weight in another run,. The conversion 

vs. (~) data for these runs are plotted in Fig. 5. The data 

r 

for both the runs fall on a single curve, and hence it can be 
concluded that external mass transfer effects are not important 
and can be neglected. The minimum flow rates employed in this 
investigation were 16 cc of hydrogen per min and 0.0143 g.mol 
of m-xylene/hr. In all subsequent experimental runs flow rates 
above this value were employed. 

4.1.2 Effect of catalyst particle size 

/ W , 

To test pore diffusion limitations, conversion vs (“jT") 
data were obtained for a finely ground catalyst powder and 
1.6 mm pellets. The results presented in Figur'^ 6 shows , 
negligible change in conversion with catalyst particle size, 
thus indicating that pore diffusion resistance is absent. 
Another criteria (40) was also applied to check the absence 
of pore diffusion. The effectiveness factor of finely powdered 
catalyst was assumed to be unity. Based on this assumption 
the effectiveness factor of 1,6 mm pellet at any level of 
conversion was calculated using the relation; 
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where ^ = effectiveness factor of the given catalyst pellet. 

r^ = reaction rate at any conversion level with given 
catalyst pellet 

= reaction rate with powdered catalyst at the same 
conversion level 

For a catalyst of given particle size, the effect of pore 
diffusion on the reaction rate can be neglected if the value of 
its effectiveness factor exceeds 0.95 (41). 

The reaction rates were calculated from the slope of 
Figure 6 at different values of m-xylene conversion. These 
rates are tabulated in Table 1 along with the effectivenss 
factor values of 1,6 mm catalyst pellet. The average value of 
effectiveness factor for 1.6 mm pellet was calculated as 0.947 
and thus pore diffusional effects can be neglected. 

In zeolites which have bimodal pore size distribution, 
two distinct diffusional processes are involved. 

(a) macropore diffusion 

(b) micropore diffusion 

Micropore diffusional effects cannot be evaluated unless the 
crystal size of the zeolite is changed. As the diameter of the 
pores in mordenite is slightly . larger than that of the hydrocarbon 
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TABLE 1 ; Comparison of reaction rates of different 
catalyst particle sizes 


Catalyst ; 1.5% Ni-HM 

Reaction ; 623 K 

temperature 


m-xylene 

conversion 

mole% 

m-xylene reaction 

rate g.mol/ (hr ) (g. catalyst ) 

Effectiveness 
factor of 1.6 mm 
pellet catalyst 


1.6 mm pellet 
catalyst r^ 

Powdered 

catalyst 

^2 

nn ^1 

^ ^2 

0.0 

0.85 

0.9 

0.944 

10.0 

0.67 

0.706 

0.949 

20.0 

0.623 

0.651 

0.957 

30.0 

0.595 

0.636 

0.937 
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molecules like xylenes, micropore diffusional resistances 
cannot be ruled out. Since the value of effectiveness factor 
is close to 0.95, micropore diffusion resistance is probably 
very small and can be neglected in the present study. 

All the sxibsequent experimental runs were conducted 
with 1.6 mm catalyst pellets. Since the mass transfer effects 
were very small, it has not been considered in rate as v/ell 
as activity studies. 

• 2 Activity and selectivity s t udies 

The activity and selectivity of nickel hydrogen mordenite 
catalysts for m-xylene isomerization were obtained under various 
operating conditions. Effect of following variables were 
studied on activity and selectivity of the catalyst. 

(a) H 2 /ni-xylene mole ratio 

(b) 'Liquid weight hourly space velocity (LWHS'V) 

(c) Reaction temperature 

(d) Nickel content of the catalyst 

(e) Nickel crystallite size 

Preliminary experiments performed with various catalysts 
revealed best activity and selectivity for 1.5% Ni-HM which 
was in accordance with liquid phase xylene isomerization reported 
by Sensarraa (34). Therefore, 1.5% Ni-HM catalyst was chosen 
to study the effect of different variables over selectivity and 
activity. 
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'^•2.1 Effect of H^/m-xylene mole ratio 

The effect of H 2 /m-xylene mole ratio on conversion^ 
isomerization activity and selectivity are plotted in Figure 7. 

The activity and selectivity decreased with increase in H 2 / 
m- xylene mole ratio in the range of 3 to 5. The activity 
dropped from a value of about 42 mole% to 34,5 inole% while 
the selectivity decreased from a value of about 75 mole% to 
61 mole?4. The conversion of m-xylene was nearly constant with 
H 2 /m-xylene mole ratio. The activity and selectivity of the 
catalyst dropped because of formation of by products at higher 
H 2 /rn-xylene mole ratio due to disproportionation reaction, 
H 2 /m-xylene mole ratio of 3 was used in all siobsequent experi- 
mental runs as this showed maximum activity and selectivity 
for the isomerization reaction. Lower values of H 2 /m-xylene 
ratio were not used as it may lead to deactivation of the 
catalyst due to coke formation, 

4,2,2 Effect of LVMSV 

The effect of LWHSV on isomerization activity and conversion 

of m-xylene to para and ortho isomers is shown in Figure 8, 

-1 . 

An increase in LWHSV from 0.5 to 4.0 hr results in a decrease 

in activity and conversion. Higher flow rates of m-xylene at 

higher LWHSV results in a decrease in residence time of the 

reactants in the catalyst bed, thereby dropping the conversion 

-1 

and activity. ' Although LWHSV value of 0,5 hr gave maximum 
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conversion and activity, it may probably be in the region of 
external mass transfer since this corresponds to an H 2 gas 
flow rate of 16 cc/min and m-xylene feed rate of 0.0143 
g. mol/hr. It has already been pointed out in section 4.1.1. 
that flow rates of H 2 gas and m-xylene above these values 
ensure absence of external mass transfer. Hence, LWHSV value 

of 1 hr (32 cc of H 2 /min and 0.0285 g.mol m-xylene/hr) was 

-1 

used. A value of 1 hr was used for subsequent studies on 
activity and selectivity because of appreciable conversion 

of para and ortho isomers. However, in kinetic studies LWHSV 

-1 

value of 4 hr was used to get differential conditions of 
conversion. 

4.2,3 Effect of reaction temperature 

In order to study the effect of reaction temperatuxe 
experimental runs were made in the temperature range of 523- 
623K at constant LWHSV of 1 hr”^ for 1.5% Ni-HM catalyst. The 
effect of temperature on conversion, isomerization activity 
and selectivity of the catalyst are shown in Figure 9. Ihe 
conversion progressively increased with temperature from 523- 
623K. The isomerization activity also increased from about 
13 mole% at 523K to about 42 mole% at 623K. The isomerization 
selectivity increased with increase in temperature from 523- 
573K and decreased with further increase in temperature upto 
€23K. The selectivity was maxiraxam at 57 3K and its value was 



Mole 
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FIG.9 EFFECT OF TEMPERATURE ON CONVERSION, 
ISOMERIZATION ACTIVITY AND SELECTIVITY 
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90 mole%. 

Plots of the conversion of m-^xylene to para and ortho 
isomers vs. reacfion temperature are shown in Figure 10, 

Figure 10 shows that the rate of p- xylene formation is more 
than o-xylene formation at a given temperature. The para 
isomer yield becomes almost constant beyond 598K. Para xylene 
selectivity was also calculated and found to increase from 50 
to 63%. in the temperature range of 523-573K. Further increase 
in temperature dropped the selectivity value to 57.5% at 623K. 

The increase in activity and conversion with increase in 
temperature is due to increase in isomerization rates.- The 
slight increase in selectivity with temperature (523-573K) 
shows that the disproportionation of xylene is not significant 
whereas at temperatures above 57 3K/ a drop in the value of 
selectivity may be due to prominence of disproportionation 
reaction resulting in the formation of toluene and trimethyl 
benzene. However, the total yield of disproportionation 
products were not more than 1,4 mole%, 

Para isomer yield was more than ortho isomer due to 
para selectivity of the catalyst. The "effective channel 

O 

appertures" of the mordenite are 8 to 9 A. The molecular 
dimensions of the three xylene isomers as found using Fisher- 

O 

Hirsch-elder- Taylor hard sphere model (41) are 7.0 A for 

, - O - ' ® 

p-xylene, 7.6 A for m-xylene and 7.6 A for o-xylene. Therefore. 



Isomer in produc 
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FIG.10 EFFECT OF TEMPERATURE ON CONVERSION 
OF m-XVLENE TO p-AND o-XYLENES 
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the p-xylene could readily diffuse out of the catalyst pores 
while the ortho and meta isomers could not do so rapidly* Thus 
the effective residence time inside the catalyst pore of the 
more bulky ortho and meta isomers would be increased and would 
produce more para isomer. A. slight decrease in p-xylene 
selectivity with rise in temperature beyond 573 K may be due. to 
increase of o-xylene dif fusivity at a rate faster than p- 
xylene dif fusivity. 

Table 2 shows the comparison of activity and p-xylene 
selectivity with other catalysts reported in literature. The 
values obtained in the present work v;ere comparable with the 
reported values. HZSM-5 catalyst gave highest para selectivity 
because of its pore dimension limitations. Ni-HM catalyst gave 
higher para selectivity compared to zeolite Y and silica alumina 
catalysts. Thus Ni^HM seems to be a good choice for isomeri- 
zation catalyst. 

4.2.4 Effect of nickel content 

The nickel content was varied from 0,5% to 5%, by weight. 

The conversion, isomerization activity and selectivity of 
catalysts with different nickel contents were studied at a 
reaction temperature of 573K and are plotted in Fig. II.- It 
has already been shown in section 4.2.3 that maxi m\am selecti- 
vity was obtained at 573K. The conversion and isomerization 
activity was found to increase with increase in nickel percentage 
from 0.5% to 1.5% by weight and reach a maxima. The values 
decreased with further increase in nickel content upto 3 wt/4. 




40 


TABLE 2 I Comparison of activity and p-xylene selectivity 
of various catalysts 


Catalyst 

Activity 
(mole %) 
and reaction 
conditions 

p-xylene selectivity 

= Tk ^ 100% 

para-i-ortho 

Reference 

Silica 

alumina 

45.4 

(755K, Vapour 
phase) 

57.1 

16 

LaY 

23.4 

(623K, Vapour 
phase ) 

43.0 

28 

HZSM-5 

7.5 

(573K, Vapour 
phase ) 

64.1 

29 


41.5 

(573K, Vapour 
phase ) 


Present. Study 


1.5% Ni-HM 


58,4 



Mole 


4i 



0 1-0 2*0 3-0 4-0 5-0 

Wt 7o Nickel in catalyst 


FIG.n EFFECT OF NICKEL CONTENT ON CONVERSION, 
ISOMERIZATION ACTIVITV AND SELECTIVITY 
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5/0 Ni-HM showed better conversion and activity when compared to 
3% Ni~HM« Blank mordenite (0% nickel) also shov^ed appreciable 
conversion and activity. The isomex'ization selectivity was 
found to increase with increase in nickel percentage from 0 to 
3% and decreased with further increase in nickel content upto 
5%. The isomerization selectivity of 1.5% Ni-HM was 90%, 

The turnover number (momber of molecules of m- xylene 

reacted per second per nickel rnetal site) va.lues were calcu- 

Icited based on the metal area calculated from XRD moasurements . 

The number of metal sites wore calculated on the basis of projec- 

2 

ted area of hydrogen molc'cule (0.120 nm per molecxile) and 
assuming linear stoichiometry of hydrogen adsorption. Hydrogen 
chemisorption studies were not conducted on these samples * 
Hydrogen chemisorption can detect the presence of small nickel 
particles present inside tlie cavities of mordenite. Since 
the size of the co.vity in mordenite is 0.8 to 0.9 nm and 
catalysts were prepared by. impregnation technique, it is 
very likely that most of the nickel particles are on the outer 
surface of mordenite. Therefore/ X-ray diffraction gives a 
fairly good picture of nickel particles dispersed on the outer 
surface. Thus the metal area calculated from XRD data is 
justified in calculating turnover niimber. Table 3 shows 
nickel crystallite size/ metal area and turnover ntomber of 
various Ni-HM catalyst samples. 

The drop in activity and conversion beyond 1,5% nickel 
may probably be due to the poor utilization of nickel particles 
present on the outer surface of mordenite. Dispersion of 
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TABLE 3 5 Nickel crystallize size, metal surface area and 
turnover number of Ni-HM catalysts 


% Ni in 
catalyst 

Ni crystallite 
size nm from 
X-ray 

Ni metal ^ 
surface 
area m^/g 

No . of metal 
sites per g. 
of catalyst 

xlO-20 

Turnover^ 
number 
xlO^ 
(sec“^ ) 

0.5 

23.86 

28.85 

2.35 

0.986 

1.0 

22.61 

29.82 

2.49 

1.136 

1.5 

19.52 

34.54 

2.88 

1.861 

2.0 

21.48 

31.39 

2.62 

1.079 

3.0 

19.18 

35.15 

2.93 

0.860 

5.0 

18.68 

36.09 

3.01 

1.282 



6 

X 10^ m^/ g 


^ Nickel 

metal surface area = 

V 



D 

V 

= Nickel crystallite diameter^ nm 


= Bulk density of nickel 
= 8,90 g/cc 


b 


Turnover number 


No, of m-xylene molecules reacted 
(metal site) (sec) 
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nickel can be affected by number of factors such as degree of 
reduction of nickel in nickel mordenite catalyst as pointed out 
by Suzuki et al. (42). 5% nickel gave higher activity and 

conversion compared to 3% probably due to high degree of nickel 
reduction leading to better dispersion. However, 1.5% Ni-HM 
gave maximum value of activity and conversion compared to all 
catalysts. XRD measurements on 1.5% nickel catalyst reveal the 
possibility of optimum crystallite size responsible for maximum 
activity. The crystallite size effect on conversion and activity 
is discussed in detail in section 4.2.5. Blank mordenite shovred 
appreciable activity because of its acid sites present which are 
primarily responsible for isomerization reaction. The metallic 
nickel formed on zeolite surface play a decisive role in the 
isomerization reaction by suppressing side reactions like dis- 
proportionation . 

4.2.5 Effect of nickel crystallite size 

The average crystallite size was calculated from XRD 
measurements and the values are shown in Table 3. It is reported 
in the literature that crystallite size may have an important 
bearing on the catalytic activity (43). The effect of crys- 
tallite size on activity is complicated by several factors 
such as metal properties, range of crystallite size and method 
of preparation. The average crystallite size of nickel particles 
have a marked effect on m-xylene conversion and turnover number 
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as shown in Figure 12. As the nickel particle size increased 
beyond 19.5 nm the conversion and turnover number values 
decreased sharply. The average nickel particle size obtained 
in the present study was in the range of 18 nm to 24 nm. 1.5% 
Ni-HM catalyst with crystallite size of 19.5 nm gave maximim 
conversion and turnover number. These nickel particles v^hich 
lie on the outer zeolite surface are so large that they behave 
like metal on other supports. Influence of metal crystallite 
size on activity of benzene hydrogenation was reported by ooughlan 
et al. (46) using a series of ruthenium zeolite catalysts. 

The effect of nickel metal area on turnover number is 

shown in Figure 13. Turnover number increased with metal area 
2 

upto 34.5 m /g and beyond this value, the turnover number dec- 

2 

reased sharply. 1.5% Ni-HM with a nickel metal area of 34.5 m /g 
showed the highest value of turnover number. 

In nickel mordenite catalysts, the formation of nickel 
crystallites is affected by the presence of acidic hydroxyl 
group as reported by Suzixki et al. (44). The growth of niokel 
crystallites on zeolite surface thus depends upon the pretreat- 
ment procedures, degree of reduction and hydroxyl groups, 

Vannice (45) has suggested that with nickel an optimum cri'-s- 
tallite size might exist, since large unsupported nickel 
crystallites have low activity. The observed variation in 
turnover number is in accordance with Davidova et al, (43) 



-xylene conversion mole 
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FIG.12 EFFECT OF NICKEL CRYSTALLITE SIZE ON 
CONVERSION TURNOVER NUMBER 



7 




48 


who indicated that there exists a definite optimum size of 
nickel particle in zeolite Y which will give maximum activity 
for toluene disproportionation. 
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4.3 Reaction kinetics and Modelling 

Though there are a substantial number of publications 
concerning the isomerization of xylene over various catalysts, 
the mechanism of the isomerization reaction is still not very 
clear. The reaction scheme for the isomerization can be 
represented as 


(i) Series reaction 

para xylene meta xylene ^ 5 = t ortho xylene 

(ii) Coupled reaction 


Para xylene 

// 

Ortho xylene • - — meta xylene 


In the first scheme no interconversion between the para and 
ortho isomers is considered, where as in the second scheme inter- 
conversion among the three isomers is considered. 


Both these reactions were successfully applied in n-umber 
of studies with different catalysts systems in both liquid and 
vapour phase. For the homogeneous catalyzed reaction using acidic 
halide catalysts, the investigation led to rather general agree- 
ment that isomerization of xylene was series reaction with no 
interconversion between para and ortho isomers (1—3). Silvestri 
and Prater (15) calculated the rate constants for the coupled 
reaction network between the three isomers for the reaction with 
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silica alumina catalyst. The relative values of these rate 
constants indicate that the direct conversion of para to ortho 
and vice versa is very small compared to the conversions of 
para to meta and ortho to meta isomers. Cortes and Corma (17) 
also ruled out the possibility of interconversion between ortho 
and para isomers for the isomerization reaction on silica 
alumina catalyst. For the same catalyst-reaction system Hanson 
and Engel (16) also successfully applied the series reaction 
model. The series reaction model was also applied by Cbllins 
et al. (28) and Chutoransky et al. (27) for large pore Y-type 
zeolite catalysts. With hydrogen mordenite catalyst Hopper 
and Shigomura (4) calculated the six rate constants accor^Sing 
to the coupled reaction scheme for liquid phase isomerization 
of o-xylene. Their results are consistent with the picture 
of a system in which isomerization easily takes place in a 
series of 1,2-methyl group shifts (i.e., meta to ortho and meta 
to para), but occurs with difficulty in the 1,3-methyl migration 
(i.e., para to ortho), Miklosy and Papp (6) found their results 
to be consistent with series reaction scheme for the isomerization 
of m- xylene using silver exchanged H-mordenite (Ag-HM). Collins 
et al. (29 ) found significant interconversion between para and 
ortho isomers. They concluded that a series reaction mechanism 
applies to homogeneous silica-alumina and large pore zeolite 
catalysts but not for the small pore zeolite ZSM-5. 
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Based on these considerations given above, in the present 
study of vapour phase isomerization of m-xylene using large pore 
(0.8— 0.9 nm) hydrogen mordenite as catalyst support the series 
reaction model would be more appropriate for kinetic data 
analysis. Since in the present work the starting isomer is 
m-xylene and there is no interconversion between ortho and para 
isomer, the series reaction can very well be represented as a 
parallel reaction. 

The parallels reaction scheme for m-xylene isomerization 
can be represented as s 

A y— R (9) 

S 

A s m-xylene 

R : o-xylcne 

S i p-xylene 

The two xylene isomerization reactions (meta to ortho and meta 
to para) occur simultaneously and independently, which made it 
possible to observe the kinetics of both reactions at the same 
time* 

Hougen— Watson formalism, based upon Langmuir Hinshelwood 
mechanism was used to derive the rate models and are given 
below for the reaction represented as m-xylene yacb p-xylene. 
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Adsorption of A controlling 


r 


s 


’"a ' 



* ^s'- Kg) + 


( 10 ) 


ainqlo site reaction controlling 

r = 


‘ 1 + K^^A KrP^ + KgPg). 


Dua l site reaction controlling 


K' (P^ - ^) 


‘l«WWKsPs>' 


(11) 


( 12 ) 


Desorption of S controlling 


(1+(K^+K*K3)P^+K3P3 ) 


( 13 ) 


who re i 


r, = rate of ■ isomerization of m~ xylene to p-xylene,- 
S 

g.mol- 

(hr) (g. catalyst) 

k* = overall rate constant for reaction of m-xylene to 


g*mol - 

p-xylene, catalyst) (atm) 
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rate constant for adsorption of A or S on 

the catalyst surface — — r—:; — rr r' "l ~~ T 

(hr) (g. catalyst) (atm) 


^‘^sorption equilibriim constants for components 

v-1 


K’ 


A, R, s respectively, (atrn)' 
homogeneous reaction equilibrium constant. 


Similar expressions were derived for the reaction m- xylene 
o-xylene. 

These rate models were simplified by deleting some of the 
terms because of their negligible contribution. Table 4 
summarizes all the models including simplified oneswhich were 
tested for their suitability. 


Kinetic runs were made at different temperatures ranging 

from 523-623K under different m-xylene partial pressures. 

-1 

Liquid weight hourly space velocity of 4.0 hr was used to 
ensure absence of diffusional resistance and also to obtain 
differential conditions of conversion. A typical kinetic data 
at temperature of 623 K are given in Appendix A. 

The effect of partial pressure of m-xylene (P^) on 
conversion and p-xylene selectivity at different reaction 
temperatures is shown in Figures 14 and 15 respectively. The 
partial pressures are the average between inlet and outlet 
values. The conversion increased at all temperatures as 
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TABLE 4 s Summary of rate models tested 


Rate ^ Model Rate equation for the 
controlling reaction meta to ortho 

step 


Rate equation for 
the reaction meta 
to para 


1. Adsorption A-1 
of ra-xylene 


^A^^a“ k ^ 


-R 

I+PrCiT 


V =’ (p 

A^ A K‘ ^ 




K 


1+P,(^ +K„).+KrPr 


S K 


2 .Surface 
reaction 
on single 
site 


SS-1 


■ -Pr 

k (p _ -ii ) 
A K ^ 


•r 


( 1+K^P^+K^Pj^+KgPg ) 




k-(P^- ^ ) 




SS-2 rR= 


k (P^-^) 


^A^A ^R^R ^ 




K' ^ 

(H-KaPa+KrPr) 


3S-3 r 


R 


^ ^^A~ K ^ 




(l+KAP^-i-KgPs') 


r = 
R 


(1 + K^P^) 




(1 + ) 


S3-4 
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TABLE 4 2 continued 


3 . Surface 
reaction 
on dual 
site 


) 


DS-1 r^= 

K. 








( 1+K P +K P +K P 
^ A R R S S 


DS-2 r„= 

K 


^^^A " K ^ 


(It K^P^ t Kj^Pj^)2 




^^A ~ K' ^ 


DS-3 r 




R 


(1 + + KgPg)^ 


^^A ~ K' ^ 

- q P 

( It K P t K P ) 
a A S S ^ 


DS-4 x:^= 


- r ) 

(1 + k ^ p ^)2 


^s= 


(1 + K^P^P 


4 .Desorp- 
tion of 
ortho/ 
para 
xylene 


D-l r. 


~ K ^ 

l + (K^+KKj^)P^+Kj^Pj^ 


^S^^^A ~ K ^ 

1 + »A+™S>WS 
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FiG.14 EFFECT OF m -XYLENE PARTiAL PRESSURE 
ON CONVERSION 



f%.a)OUi) A^iAipai as auajAx-d 


I- iG-15 EFFECT OF m-XYLENE PARTIAL PRESSURE ON p-XYLENE 
SELECTIVITY 
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the value of increased whereas the para selectivy remained 
almost constant. At higher values of P , the concentration of 
m-xylene in the feed stream is high and hence the conversion is 
high. As the value of P^ increases the relative proportions 
of ortho and para isomers also increased. Since para selec- 
tivity is the ratio of para to (Para+ortho) and the reactions 
(meta to ortho and meta to para) are similar in nature, 
the para selectivity remained constant. 

4*3.1 Mathematical modelling and data analysis 

A non-linear regression algorithm utilising Newton-Gauss 
technique was used to obtain, a mathematical fit for various rate 
models. The program minimized the residual sum of square (R.S.S.) 
during the regression ^ The nonlinear computer program basically 
improved upon the initial estimates of various constants of 
the rate equation until the R.S.S. could no longer be reduced. 

The starting values of the parameters were estimated by linear 
regression. A computer program listing is given in Appendix B. 

The model discrimination was done based upon the requi- 
rement that kinetic and adsorption constants had to be positive. 
Some models were rejected on the basis of improper trends of 
these constants with temperature. The data were most satis- 
factorily correlatedby model DS— 4 (surface reaction on dual 
site rate controlling). Although/ this model has close values 
of RSS compared with SS-4 (surface reaction single site rate,- 
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conti'olling ) , the choice was based upon the requirement of 
proper trend of the constants with temperature. Also, the 
average absolute percentage error for model DS-4 was lower 
than model SS-4. 

Experimental rates were plotted against theoretical 
rates calculated from model DS-4 for para and ortho isomer 
formation in Figures 16 and 17 respectively. The experimental 
rates agreed well with the predicted rates as clear from the 
Figures. Table 5 gives the values of rate constants and 
adsorption equilibrium constants of the model DS-4 along with 
average absolute percentage error and R,S«S. for different 
temperatures , 

The effect of temperature on the rate constants for the 
isomerization reaction is shown in Figure 18 in the form of 
Arrhenius plots where the rate constants are plotted against 
the reciprocal of reaction temperature on a seiini log scale. 

The activation energy of the reaction was calculated from the 
slope of the straight line. The values of activation energy 
for para and ortho isomers from meta isomer were found to be 
nearly same as expected because of the similar nature of the 
two reactions. 

Comparison of activation energy values obtained in the 
present study with the reported values is shown in Table 6, 

The activation energies obtained in the present study was 
slightly lower than the reported values probably due to micropore 
diffusion limitations in mordenite support. 
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fig-18 EFFECT OF REACTION TEMPERATURE ON 

REACTION RATE CONSTANTS (Arrhenius plot) 



TAi'L.i 6 2 Comparison of activation energies with reported 
Vctlues for the isomerization reaction 


Catalyst 

Activation 
energy kcal/mol 

Reference 

Hydrogen 

bromide/.Muminium 

Bromide 

22+ 1 

2 

Silica alumina 

25,5 

16 ■ 

Hydrogen 

morderiite 

27+2 

5 


34-37 

6 

lii-Ml’i 

21-23 

34 

i: I X *■’*» i' .. 

15--18 

Present 


study 
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CHAPTER 5 

CONCLUSIONS AMD RECOMMENDATIONS 


6.1 Conclusions 

Nickel supported on hydrogen mordenite was found to be 
an active catalyst for vapour phase isomerization of m-xylene. 

Conversion of m-xylene remained almost constant with H^/ 

m~ xylene mole ratio (3 to 5), decreased with increase in LWHSV 
-1 

(0.5 to 4 hr ) and increased with increase in temperature 
(523 to 623K). The isomerization activity decreased with E.^/ 
m-xylene mole ratio^ LWHSV and increased progressively with 
temperature. The selectivity decreased with H 2 /m-xylene mole 
ratio and showed a maxima at a temperature of 573K, 

Increase in nickel content of mordenite (0.5 to 1.5%) 
increased the conversion and activity and showed a maxima 
around 1*5% nickel. Blank mordenite gave better activity and 
conversion compared to catalysts having nickel content greater 
than 1.5%. 

A pronounced effect of nickel crystallite size (18-24 nm) 
was observed on m-xylene conversion. 1.5% Ni-HM catalyst 
having crystallite size of 19.5 nm gave maximum value of 
turnover niintoer. 

The kinetic data were most satisfactorily correlated by 
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represented as 
Pq 

_ - ic? ) 

~ ^5 (meta to para) 

k(p^ - ^ ) 

- 7r (meta to ortho) 

(1 + 

The activation energy for the isomerization reaction was 
calculated as 15-18 kcal/mol. 

It is expected that the present results will contribute 
to the understanding of kinetics and the role of nickel in 
Ni-HM catalystfor m-xylene isomerization, 

6,2 Recommendations 

1. Ciiaracterization and morphology of nickel particles in 
Ni~HM catalysts should be studied using other techniques 
such as Hydrogen gas chemisorption and electron microscopy. 

2. A conparative study of Ni-HM catalysts prepared by ion 
exchange and impregnated techniques should be attempted to 
the role of small nickel particles entrapped inside the 
mordenite cavities* 

3. Effect of catalyst preparation variables like calc'ina'tlon 
temperat\are/ reduction temperature & time should be 
studied to establish the nickel crystallite size affect 
on activity and selectivity. 


the model 


and 


■R 
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4 . Deactivation studies should be carried out with nickel 
mordenite catalyst to check the decay of activity with 
time . 
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APPENDIX B 


COMPUTER PROGRAI'l FOR CALCULATION OF RATE AND 
ADSORPTION EQUILIBRIUM CONSTANTS USING LINEAR 


AND NON-LINEAR REGRESSION ANALYSIS 
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•1 


YCK)=:l 

XtKrU 

XtK,2j 

X(K,3) 

X C K , 4 ) 

XtK,5) 

XtK,6) 

Ns 4 

(lU TO 

y(K)=i 

XtK,U 

X( K,?) 

X( K,3) 

XtK, A) 

XtK,5) 

X(K,6) 

ds4 

no TO 
YCKJ = 1, 
XfK,l) 
XCK,2i 


.0/HP(K) 

= .l .P/ABtK) 
=FSCK)/AB(K) 

= PSlK)/AmK) 
SPR(K)/APCK) 

=0 . u 
=0,0 

16 

.0/HP(K3 
= .l.U/ARCK3 
=PA(K)/AP(Kj 
=PR(K)/APCK) 
=PS(K)/AR(K) 

= 0,0 

=o.u 

16 

fj/(RRCb)*=«'0,5) 

1 .0/CAIi(K)t*0.5) , 
PACK) /(ArifiO =**0.5) 



■'lb 


4 1'l 


5t> 


6t» 


7 b 


XtK,3) 
X ( K , 4 ) 

X t K , 6 ) 

Ms 4 

<;i( rn 

Y(K)si 

xtK,n 

XtK,2j 

X ( K , .1 j 

X I X , 4 } 

•n S 3 
'10 TO 

XIK,1 j 
X(K%2J 
XCK,2) 
X C K , 4 J 
X ( K , 5 J 
XIK,6J 
M = 'l 
OU i’ri 
YtK) = :t 
XCK,i j 


-P,S CK ) / c AO f K ■) ■■*•♦0, 

“o.u 

= 0,0 

1 h 

.0/KRtKj 
= 't ,0 /hH{KJ 
=:PUK)/An(F) 

= PACK)//0HK) 

= P''HKJ/aR(K) 

— ( _ ii 


1 f) 

. i‘ t K > 

= I .O/iOHK) 

= P.'1CK j/ABtK) 
=r'R(X)/AA(K) 
=>},()■ 

= b , 0 
TM , 1) 

1 b 

.O/RIHK) 

= 1 ,0//ib CK J 


XCK,2 J=:bSCK)/AB(K j' 
XtK |)=pntK)/ARC^j 

X(K,Bj=:u.O 

XCK,H;r<,,,; 
rjs i 

Gl-! TO 1 0 

YtK)=l.()/HR(K) 
X^4) = i .U/A«tK) 
JtK.2)=!PS(KJ/ABCX) 
X{K , 2 ) s(j y 

X(K,4)=.,).y 

XCK,5J=0.0 

XCK,6)=0.0 


Ws? 

flu TO 

Y(K)s:l 

X(K,1} 

X(K,2) 

XCK,3) 

X(K,4) 

X C K , 5 ) 

XCK,6) 

0 = 3 

Gt) TO 

YiKJ=J 

XtK, i ) 

Xl»%2J 

XCX,3) 

XtK,4 J 
XCK,S) 


J6 

,0/RH(K) 

=i.o/abckj 

=PA(KJ /ABtK} 
= FR(XJ/AB{I<) 
= 0.0 
= 0.0 
= 0.0 

16 

.0/RRtK} 
=1,0/mB(K} 
=PA (K)/AB(K} 
= PvSCK)/ARCX '3 
= 0.0 
=o.n 





l'j> 


9 b 


lU6 


3 16 


126 


!b 

3.2 

C 

200 

RiM) 

R 0 1 
C 


XlK,n)::A,,i 
S i 


Ou TO 3b 

XCK,U=;3 .u/aB(K) 

yCK,2J=(:'AtKj/A9{.Kj 

\(K f ^ ) si: i) 

'CCK.'lJr.uO 

>’ t K , S ) =b . 0 

X(.r,AJ=:u.O 


Ob ni 1,6 

Y1KJ = .1. .0/CI^3<(6)4*0,b3 
XIK,1 J = l.u/^.Abf^3+»^),5) 

X ( K , 2 j = p A ( K ) / ( A I j f K 1 ♦ ♦ 0 , 6 ) 

X t K , 3 J = P R C K ) / ( A b ( K 3 * * 0 . 5 ) 

X(K,43=0,0 

X(K,Si=i),0 

XtK,6)=0.0 

H = 3 


GO I’O 3. b 

YtK) = l.y/(;RRCh )=^*0.53 
X(K,l) = 1.0/(A!i(K')^*0 -5) 

X (. K , 2 3=:f»ACK J / i Ab ( K) 5 ) 

X t K , 3 3 -PS t K 3 / ( AB ( K) 5) 

X(K,4)=u,0 

XtK,5)=0.0 

X(.K,6)=<},0 

M = 3 


GO TO 36 

X(!<) = 1.0/CRR(K31‘’t'0,53 
XCK,3 ) = 1 .0/CABfK)**t).53 
X ( K , 2 3 sp A ( K ) / ( Ab f K 3 ♦ ♦O . 5 3 


XCK,33=0,U 
XCK,4)=!0.0 
XlK,SjsO.O 
XCK, 63=0.0 
M=2 


GO TO 16 
Y(;K3 = 1.0/(iHCK) 

XCK,U = 1.0/AB1F3 

X(K,2)=PA(K)/AR(K) 

X(K,33-FA€KJ/AP.CK) 

XCK, 43=0.0 
XCK, 63=0.0 
XCK,6)=b.O 
M=3 

G(j TO l6 

coNTrsint; 

CUMTIWOE 

rOLu SHARMACY,X,PA,PR,PS, 
CkhU SHARrtA(Y,X,f^,W) 
9M=NN+1 

TF(Nrt.GT.143 GO TO 200 

GU TO 400 

CO-m-MllE 



6T0P 

EU!3'')ATf /tOX, 'PRUBGFM wnUBEP 
Fi.ROATC/lOX, 'MOUEG NOiffiER 
Rpn 



T2/3 

12/3 




V 

r 


c 

1 0 1 
1 1 • 4 

1 


IB 

r 

r 

5 

6 
C 


7 

C 

1u 

r 

11 


u 


14 

C 

c 


15 


B2R 

Rio 

R20 


roj^KOUTiMli SHarImA C Y, X , P/\ , PR , P,5 , ’i , 

n ^ v*' f '^1' 1 » X f 2C , ?i.- 1 , X A ( 2<S ?•. 1 

1 X T i ( ■] I J , Y Y ( 5 5 ^ ^ ^ 2 <s 1 , vH 0 ER C 2 0 1 , X 1’ ( 2 ('■ , 2 u 1 , Y i Y i ? s * , 

PT.‘Uv4oTu''i R C20i 20 J ,f]1 t2U ,201 ,aFRR( 2«: 1 
KiOE.AT OF OUTPin- STATEOKf’Tr, . », v.. i 

r,o"iiV’T-lt '*7*-’^'''^ *' VECfOn T^'J 
YV( li=:Y(TJ 
nu IR I=:J.M 

no IP j=i,o 

XA(I,01=X(T,,1j 
CEr,cnbATTUN Of iW.iVfs 
nEl'Et''iTnE THF TKA.IRpncjR OF X 

ini P 1 = 1,0 

I'tJ b J = 

XTCO,Il=XtT,.Ij 

Fi;Fi,.> XTXf.fJX[0=XT(NXw)*X(MXM) , 

Ot ( / .l = l,w 
Ot.l / 

XTXtl, J)=0.(.5 
IHl 7 E=X,f.| 

Xl’Xa,>T)sXTX(i,U)+X7f i,K)5t:X(K,.T) 

TtJ FI,jn TNVERRF OF X^X, 

CAI,L -IATToV CXTX,i'!.b1 ,0 ,r)ETfc;R0V 
FI’jO X.'^Yf '•OsXTf bXi-' W-Y (.mj 

'■'0 11 r=i , 0 
XTY(T)=:O,0 
011 121 = 1,0 
nu 12 u=i,M 

XI Y ( T ) SXTY ( 1 1 +XT ( I . J ) IJ 3 
FIOD b(lO=XTX(MXN)XxTY(Al 
00 U 1=1, M 
1HI)=0,0 
Oil 14 1 = 1,E' 

no 14 js-i ,n 

iUTl=i.d)fXTXCT,.1)*XXY(J) 

PtJ f j; ) =B ( T ) 
rOMTTi^UE 

PKTKT82U, CB(I) ,1=1,0) 

Bl N1'P=Y'.1FAN-BC 1 J *X1 HF-R (2) *X2(V|M-.R(. X )*X .X -iM- 
lHC4)i1X4MM-af5)*X50i'j-Bf63*XoY,v 

SSR= 0 ,n 
AHFUSO.O 
nu 15 1=1, M 

yi!AT(I)=nCl jl'XXCT/l )+uf2)*XX(l ,2)+RCi)*XX(T,Xj + 
lK(4)*XX(I,4)+B(53*XXlI,5)+Bfb5*XXf i.,R) 

EBHURCIl = Y(l)-YH»'f(l) . . » 

ABFk=AnER + fERHnHCI)/Y(T)) 1^100,0 

SyFH(i)=FHRUR(T)’^*2 

.SSR=SSRt.SgFUCI) 

conttool 

AVAaSE=AREP/FEnATiMj 

PHifiT R25,SSR 

PRTMT Bia,AVAl,4SE , , 

FUPi4AT(/1 OX, "BESTDUaT. SUM OF SOUARF 
FOMMATClOX, *AVFHAGF AiiSur.UTE ERHOh 


,E1 2.5/i 
,E12,5/J 


FuRMATC/lOX, 'FINAL PARAi'iFTRR VALUES ' //I OX, 0 ( iX , FI 2 , S I) 


2C J , 





n 

r 


C 

1 G 
< 5 

'’G 

41; 

45 

5 (i 

7f: 

«0 

«f> 

-’t; 

05 
1 C' ‘ » 
i(;5 

1 !'• 
r 

til 

14I' 

t5''’ 

IbH 

1 70 
70<> 
20*^ 
2iO 
22'' 
2i0 
2b (• 
260 
270 
C 

^10 
4 2C 
BiO 
3 40 

350 

355 

BoU 

370 

C 

3 MO 
3'iO 
4CO 
420 
430 
450 

455 

460 


OETlJiG; 

0:if) 

5 H 1 ■ fft n ^ lih „ "'I; ^ f A ' ' DFTKRM j 

■601 I Ink for F:A’^HTX IWreRRinN 4'ID .>ni,nj.T-iM hf 
ri"CiO';FfUO‘iT O/'^trTX UF SIZF. w 
vfctup flF ST'aF '•! 

' ■•MiFFOlK UF VRCTi.lRR 


iOLTA^'FnvjS r,,'G» 


IF 


'1 I 


;b'T Ti,i ZFO.n,f!.'iLi IWVFRRt IS CU-’OH i Fu 


O i OTMOijfUi?. 


i»hTEtM-»vA6HF UF Oh/l'i-.’I’.vtTia Ai'iT FETURmFu iW TnF '6,0'T!;iF 
0 i 0 A ’••( i? T I.i 0 A ( 2 0 , 2 0 ) , 0. C 2 0 , 2 0 1 , 1 0 T v 0 T ( 15 0 > , I '‘i 1 ) F X ( il , 4 j 
F2ni"ALt4CF ( tinjW,.TK'n„ Ft , f ICiO.OM^.ICniOOi) , f 

T.;n;TrJ,,i7*rr:[nH . - 

I'h'^FO 1=1 ,<) 

00 20 J= 1,0 

T0H/f)t(J)=;O 

FEAkro FHk PXvnT h'LEOEWl' 

Dll !i50 (.si,. >■4 

Ai'!AX = fi,0 

DU 105 J=l,,a' 

TF( iPXVfiTCJ)-!) 60,105,o0 
Oil iOu K=;1,D 

TF( rpcVi,i'!’CK)-n80,ton,7.*o 

r F r A 1 A V - n H 5 t A ( J , F j 3 ) 0 i , l 0 m , 1 0 0 

ikn.i=.i 

rCf'10M=K 

A.-iAXsABSCnCo ,(\3 ) 

CIOMT.OIF 

CtfOrT'jMfc; 

TFT V n . 1 - ( JL c U TOi M) - 1 p i V ( J T a C n jj U . 3 ) + 1 
T,.,TKt'C!lA’'jGR PtfWS 'in PUT PIVuT FyFkF ;T 
T t ( 1, F' U . ’ -ICO 1,1 f ! X ) 1 4 ( 1 , 200,1 4 0 
DFTKO ir-nKThRj., 

Du 2 if ,fi™ 1 , '■! 

.5riAP=Af ,IPUV.I,T0 
ACTRn.i,i,) = A(TCOLlltv>,L) 

A(IC'''1,h:‘,,T43sSwM' 

IF (I'H 266, 260, 210 
nu 25 U L=l,fi 
.5i’JAF=6ClPi)W,I0 
fUTKn'.»,jL,)=lHTCnL!if!,L3 
0( TCniO'-i,t.3=6h-AP 
TkDK,XtT,13 = lPUW 

Trin6XCT,2) = iruT,uo ^ ^ ^ 

" i V 1 1 ' F p T V n T h . n m BY p p; o t e i, k a e u t 
F lVuTsA CTCnuOi'i, ICUT.UM) 

ni'TFR.lsDFTFHf^ + PIVOT 
A(xcn6ni.i,icuLuM) = i.u 
nu 35o b=l » R . m 

AtTCnbi!tpL)=:AlTC0bnM,L3/PiVuT 

TF(.M)3l?0^iB0,360 

f^rrcobnp)7b5=PCTcnLno,L) /PIVOT 
Rtnorf.; MnwPTVfn' Hnu’.S 
no 550 ul=1fp 
TKfLt-rcnLDMl 4 00,550 , 400 
TsAlIil , .LCULlfMj 

AlTil.inUTiUMjrO.U 

D(.l 450 , .m 

AIt4l,L) = A({a,b3-AlTCnbnH,L)FT 

- TF(m5550,550,460 
no 500 b=l ,4 



S|tn 

55 f' 

r 

60 ') 
6 i'- 
620 
6 AO 
GA<‘ 
6 b',’ 
6 o" 
fi 7 f> 
70 " 
70 S 
710 


1.1 

1 2 

001 
7 t':'’ 


P 1 ^ ) “0 t TGOLHo , t, ) 4^? 

TuTt^RC^oGC CnLU/H?«S 

S !'i ^ J RO W ) 

t*< ,..K:OlOOn:;SwAO 
Ch'i-n i’lt-, 

CU''’Tr,.("K 
"li 1.1 KsT.n 


n-( LinVirCK j ,(:,p 
ruMPxr,r»f., 

OKTijn,^ 

POI (T 091 


Jj an T.J 12 


103 oTwOuLAP'/J 

‘< 1 ,. I UH 



c 

c 


7U 

A 


122 

121 


r» 

r* 


Non Linear ReGiRessjom 
Anav-ysis 

i>i vii';;us.r,i'i Av2,2j 
r'ATi;.Pi4AL, fjSgP'i'iM.T.SOi^nN 
Y,T,i%OC 

fii.»ro4M+ . iti'j (T-’-l )/? 

Pg r.K=!l,5 

'■’' = 2 

.'g,AiJ’«!,p;gC 

<^7S^ 

l.Y=2 
LlM-l 
1 1 *'( — 1 

TVPi'jf.RjC 

Pit 20 Isj.M 

ai'.AU*,Y(.T),(TtT,Jj,J = l ,3) 

COHI’TOUK 
011 3 1 = 1,11 

11 J,CTtl,JJ,j = i,3) 

IiOU01 = l 

'■''hXCAL=4O0^.<4 

[•:ta=o,5 

X 1014=1 f,0*SgRTCX02AAF(Xi‘nt,)) 

r.TFpMX=10.0 

PHI>IT 7 5 

P.vIOT 85,tXtT),l = i,N) 

CFAI0=1 

C. Ah h EQ 4FCF ( M . N . ISOl'nw , OSgMuN , I PRi M 1 , rt A XC A h , S 1 A , X ‘ft t .j , 
1 X , KS tj ItiO , F y RC , F J AC , J , 0 , V , 0 V , W I TER , F , 1 vi , I,i i F 

PhloT lul ,TFAlf, 

PKIhT 103, (XCu) ,J=1 ,H) 

CALL LSvGHDtM,W,FVFC,FlAC,LJrG) 

DU 121 1=1,0 
DO 121 vJ=l,0 
All, J1=0,U 

Dll 122 K=1,M , , 

A ( T , Jl = A( X ,U ) +FJAC (K, 1 ) ^FjAC Ctl ,U ) 

CoOllriUt 
PRTOT 120 
RO.S=t',D 

soy=o.o 

no 40 1=1, M 

Uc.iM=(l.0+X(2)»T(I,l J) i ^ , 

AOU = X ( U *X C 2 ) ^ C t C 1 , 1 ) -T C 1 , 2 1/FOC ) 

Dr:M= C I . O + X ( 2 ) *T q , I J ) * ^ U 1 ^ ^ 

AOU = XX I i *X C 2 )^ ( 1 C 1 , 1 ) “T X X q j /FOC ) 

nt,M= X 1 ,0+ XX C2 j /t:QC 5 »t ( i ,2 1 ) i ( 37 j 

AivO=XCil*XTCI,l)-l'Xl,2)/KyC) UM/) 

S5=A0U/DFM 
FVFC(i)=SG"X(I) 

.Sli()=S5Qf ( FVFC X I ) 2SS ) ♦! OU 
RiJ,S=Ri>S+FVEC(X)*'»'2 

01^10X125,138, YX I },FVEC Cl j 


.1 i ) 


,:>TfcP,..X, 

ATi,) 


i (MlO) 
1 (tMl ll 





f) l» 

1 U 1 
1 0 ? 
i if ? 

I li -i 
lijS 
i if(i 

II '> 

I 2 f- 

I 


Ci.t!! t’T.iil,-; 

I ij r 1i)2,RSS 
‘’iUrtT 

rrHi ^ ^ ' '' ' ^ 

K.IR !<VC(10X, 'Ii<T iTaT. 

? . . R .1 A I ( i lU , 2 f 2 X , f-' 1 2 . 5 J ) 

’• 1, i n 'i T ( / / K' A , * F A H IJ R E a I i 
K.,,,,MA .,'(// 1 (‘.X , ' Hu t.Xit Sg Hi- 
>*;n-igA,i, (//iOX, 'A‘i' XHt; Fnx?-^T ' 
' ‘ :’C//1'U, 'T.iF FSTx'IaTlD 


CliEoSS r’nix 
TYRF ' 


HEvjI' i'.So T u '•’ * / J 


' r xi// f 
isOUAuFb' 
,21x15, |>//J 


iA . 


nKAuIi,.?' I 
A he '//j 


1 




(V/IOa, 'TaE 
F>.<'A,iA i'(/Kl5 ,B J 

FuRilA rC/luX, 'AnijOLUl'E AVEAAb 

l;'li4uyi(SX,Aij( ' “IOjEu VAuHix* ,5 a, 

1 5 X , iJ X r Ei’.iilxMC F ' / 5 X , 6 u ( l M » ) j 
Fi*RMM'C}fU,4(4X,Fo.D) 1 

!; .* f * ' 1 


4 ti 1 g , i i / / ) 


EHRijR = ',El2.5/) 

71 L TT>,. # E! t^ 


sP'X. 


i.,. 

r» 

V:" 

r 

c 


2 t< 


c 

c 




2 U 

40 


4,ji\lUluT111.2 LSyFufl 1 IFLaG , '1 ,>! , XU , b ¥,;,rC , A 4 , TilX , , f. • ) 

Oir.r 'ooriiii FVHCCt20),Wt4i),XC(2),;i^vfl) 

DWiv'loIiiM 7120) 

CuFi'tOi.i !t,i',KOC 

Du 2A< 1 = 1,''. 

,')i.i':u'i=(i,o+xri2jii''iri,i) j u^no; 

AMUi-)=XC c U *xr (. ? J * XT ( 1 , 1. J -f C 1 , 2 ) /E^C j * ('-110 ) 
Di 2(=(1.0+XC(2)4T( 1,1) j*'??. iCilU) 

A ; tt| = XC C I ) «XC ( 2 ) ♦ ( T C i , i ) -T C I , 3 } /fcOC ) I ( Hi i ) 
r'ow=(i,o + cxc( 2 )/Eyr)*TCl, 2 ) j Um 7 ) 

.AOnsXCCl )4C I CX,1)-TX 1,2 J/EQC) lC-i7) 

2os=Aon,j/oKj,nM 

F7ECCCI>=i<E.’>~7Ca) 

co4xr.jHh 

uhTuno 

K,(D 


OOUKnijT I Ot, ijSyMOH t M , M , XC , F ECCC , F JrtCC , OOU , S , XGhf. uc , 7 1 Ti2P , « F , 1 • 

I L I W , vf , ti i4 J . . 

niMLM.'ilUH FJACCC20,2) ,F7F:CC 1 20 ) , S(2 ) , C 41) , XC i 2 j , i ’.■• t ,1 / ,G(5C'i 
Flil{,4f;y=F0lDEFiFVFCC,FVFCC,M) ^ 

CAOL OSgGKDCM,N,FVECC,FJACC,L.jC,G) 

GTGsFulOFF (G,G,N) 

R E T’ U R H 

Kf.n 

sIm}^unFirLSVGRDC^7j^7F7ECC7FuAcC7nuC7G' 

niM£,NoIUO FJACC120,2) ,FtfF.CC{.2U) ,G<.2J 
f>U 40 v)=1 ,M 
5UM=0,0 

DU 20 T=1,M , • 

SMM=f>UM + FJACCtI, J)*FVECC(X) 

CUMTIrtOh 
G(,3) = 6nu + fill!4 
CUNXlHUE 
RFTURi^ 

RivD 



^ ^ ^ ^ r ^ J ) I f’F 1' ? 5 , s ) ; 


F i: ( /i T 


1 ij 


2.J 


nj2 
1 b 


8 

9 


11 

13 

14 


17 


S;.ti'’)i<nu f'iai:; 

r.'Dtc=-.i 
'^’ii io r=i ,9 

•>u iO = 

I I i t t ,.Tj=!FFt 1 ,,j j 

•) , ■' 'VI 1 ^ ^ ^ nu Tc , 5 ) 

• > i (i«*i 1 ^kGov)/ ( t Lu A i' C M" W ) J 
in, lb t=i,9 
'■)»/ V> 9 

Tf (A,t'.;9,jj un TO 20 

I’.i! rn 1!) 

Fj TCU>f)s:GHiMAl%FK7(T,Jj 

.3l.;7=KFTtl ,JJ 

but = SviH l'( ABi>(FF'l(i,j) J) 

3o:3,C,XCTJAGl)l,KGbO 

!• I .H. uVl t i<>A, * X t ' , I i , ' ) = ' ,&1 b , a , 3X , ' GDIS ' , £ t 5 . o , iX , % £1 D . o ) 

C,i ii'J r T H'lF 

«l,Tu9,^ 

GmI) 

GuriHOuTinfe HAlIiiVtil, A,X,T9D1C,NHC) 

:UMc:'‘lGigfi TROsKb) fJCOLCb) ,JUFU(5) ,Yi5j , Af.-'iRC.viPC) ,^(.'1 j 

-iBSs i ,0F>15 

•■•Ax=n 

It (w.aE. 2'):) Gii TO 5 
GO TO lb 

0£T£R=1, 

no 18 K= 1 ,N 

KFlIsK"! 

P1VUT=0.0 
no 11 Is:i,N 
no ,11 j=i.N 
IF CK, £0.1 I GO TU 9 
DU <i iSCAi'i = l,KMl 
OU 8 JSCAiX = i,KMA 

TFCl.ftO.Thnu CIRC aim)) go Til It 
IFCd.oO.JCOLCjSCAiM)) GO Tu il 
COiiTlMUF 

TFCaBGCA(I,J)) .LE.ABS tPIVOT) ) GU TO U 
P1VUT=A(.T.J) 

TaQ*»(K) = T 
JCnG(K3=J 

CuMTli'aiO 

IFCARSfPlYn'D.Gl.tiPS) Gu TO U 
GO TO 7b 
IHOi^KsIROwCK) 

JCOGKsJCOLCK) 

DFTFRsnFTLP+PjVUT 
nu 14 J=1,MAX 

At IR0i(vK,vI 5 “AC IRuWK.J) /P i VuT 

rOMTTNfIK 

..UMII JC0S,,K) = 1. /PIVOT 

nu 18 1=1 

AIJCK = ACI,OCnGK)^ 

TFCI.FO.IiirojK j GO TO 18 
ACi;JCnuK)=“Ai.TCK/PlVOT 

XFCjlNElk’OGK) ACl.dl-ACI^Jl-ATOCK^AClPu'-FGOl 

cgntibue 





7. 1 

2d 

27 

2d 


29 

75 

45 

301 

302 


no 20 Isl ,^' 
iHi:ua=iROA( n 
jcoi,isjcnu( i 3 
JORtXXROfcDsJCODT 

If ( IMDIC .GB; , 0 ) A ( jcuf I ) ~ a t THOtt T , max j 

I.WTCH=0 

Du 22 Isl-Mi-il. 
tr;i=t + j 

t)U 22 0=:IP1 »i 

It CuruCiCOJ .GE.JuROf i ) J GO TU 22 
Jj;EUPsJuRu( J) 

JuRi.)(j)=auPD( J.) 

JuR'iif DsUTF'rtP 
Tin'CHs'IivTCHtl 

1 F { 1 Fri' r H / 2 * 2 . i j R . 1 0 T ( : H 3 i ) f r E l , - - n & T F R 
IF( jL^iuIC.iiE.O) GO Tu 26 
GU TO 75 
00 2« J=1,iV» 

Du 27 1=1, D 
IROif»J = IRnwf i) 

JCOi..I=JL0i,(l) 

YCJCOoI J=A(iRuWj.,j) 

CuMrirtijft 


C 1.1 M '1‘ 1 i 'i ( J F 


Du 2H T = 1,.i' 
ACI,J) = i ( i'i 


DU iO 1=1 
DU 29 J=1>P 
TRnwJ=IROv« (u) 
JCObJs:JCOL(J) 
YaHCU'iJ3 = A(i,JC0ijJ) 

CU^TlwUc, 


DU 30 J=1,H 

ACI,J) = HU) 

PR I iff 301 


DU 45 1=1 ,n 

Print 302£(ACl,J)f-I=l,N3 

FURMAT(/10X, 'UAl'RiX TmVFRSF IS CAuDuD /J 
FURMATCI(.'Ari{4X,El5,83/.) 

PUTURu 

END 



CHE - iqSH' ri-5R-G-yAP 



